Pattern Formation in Exciton System near Quantum Degeneracy 
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We discuss models of the modulational instability in a cold exciton system in coupled quantum 
wells. One mechanism involves exciton formation in a photoexcited electron-hole system in the 
presence of stimulated binding processes which build up near exciton degeneracy. It is shown that 
such processes may give rise to Turing instability leading to a spatially modulated state. The 
structure and symmetry of resulting patterns depend on dimensionality and symmetry. In the 
spatially uniform 2d electron-hole system, the instability leads to a triangular lattice pattern while, 
at an electron-hole interface, a periodic Id pattern develops. Wavelength selection mechanism is 
analyzed, revealing that the transition is abrupt (type I) for the uniform 2d system, and continuous 
(type II) for the electron-hole interface. Another mechanism that could possibly drive the instability 
involves long-range attraction of the excitons. We illustrate how such an interaction can result from 
plasmon wind, derive stability criterion, and discuss likelihood of such a scenario. 

PACS numbers: 



I. INTRODUCTION 

Recent experiments observed striking spatial pattern 
formation in two-dimensional cold exciton systems in 
photoexcited AlGaAs/GaAs quantum well (QW) struc- 
tures 4 ^^ 5 .. The photoluminescence (PL) patterns span 
macroscopic scales in excess of 100 /im and include con- 
centric ringsi* 2 * 3 ^*^ as well as 'bright spots 'ii In addi- 
tion to that, the electron- hole system exhibits an abrupt 
transition at ca. 2 K in which the outermost ring 'frag- 
ments' into regularly spaced beads of high PL inten- 
sity^. While the gross features of PL have been ex- 
plained within the framework of a classical transport the- 
ory, attributing the internal rings to nonradiative exci- 
ton transport and coolingi and the outermost rings and 
'bright spots' to macroscopic charge separatio n 3 ! 4 ! 5 , the 
origin of the instability remains unidentified. 

Spatially modulated exciton density is not to be ex- 
pected in QW system designed so that excitons inter- 
act repulsively as electric dipolesi and thus do not form 
droplets 6 . The macroscopic character of ring fragments, 
of 10 — 30 fj,m in size and containing about 10 4 exci- 
tons each, as well as the observed spatial modulation 
of the exciton density which appears abruptly at a cer- 
tain temperature, calls for an explanation involving a 
symmetry-breaking instability of a homogeneous state 
to a patterned state. Such behavior is reminiscent of 
the instability predicted by Alan Turingi to occur in 
a reaction-diffusion system. Such an instability of a 
reaction-diffusion system, first predicted in 1952 by Alan 
Turingi, is characterized by an intrinsic wavelength re- 
sulting from an interplay of reaction and diffusion pro- 
cesses. The Turing instability is known to occur in cer- 
tain chemical reactions&SiiS, as well as in biological sys- 



tem; 
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In this work we propose a mechanism, based on the 
kinetics of exciton formation from optically excited elec- 



trons and holes, that can lead to a Turing instability 
in the exciton system. Interestingly, these kinetic effects 
build up and become especially strong in the regime near 
exciton quantum degeneracy, due to stimulated enhance- 
ment of the electron-hole binding rate. The transition to 
a state with a spatially modulated exciton density reveals 
itself in the spatial PL pattern, and presents a directly 
observable signature of degeneracy. 

Although, in itself, an observation of an instability does 
not constitute unambiguous evidence for degeneracy, it 
may complement other manifestations discussed in the 
literature, such as changes in exciton recombination 12 
and scattering^ 3 , rates, in the PL spectrum^, absorp- 
tion^ and PL angular distribution^. While linking the 
observed instability with degeneracy is premature, our 
main aim here is to present the Turing instability from a 
broader viewpoint, as a novel and quite general effect of 
quantum kinetics that can help to identify the regime of 
Bose-Einstein condensation (BEC). 

Besides the kinetic mechanism, we consider another 
mechanism based on the assumption that the interaction 
between excitons is of an attractive character at large 
distances. The origin of such an attractive contribu- 
tion may be in the force on excitons due to plasmons 
or phonons generated at exciton recombination. We de- 
scribe a simple model for the plasmon wind effect, an- 
alyze the stability criterion, and consider the possibility 
that the observed periodic patterns are related with such 
an attraction. 



II. KINETIC INSTABILITY 
A. Model 

Following Ref£, here we consider a transport theory 3 , 
formulated in terms of electron, hole, and exciton densi- 
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ties n e ,h,x obeying a system of coupled nonlinear diffusion 
equations: 



(e) d t n c = D c V 2 n c - wn c n h + J c 
(h) d t n h = D h V 2 n h - wn c n h + J h 
(x) d t n x — D x SI 2 n x + wn^n-^ — 7n x . 



(1) 



The nonlinear couplings account for exciton formation 
from free electron and hole binding. Here, no attempt 
has been made to describe the detailed and complicated 
density and temperature dependence of the physical pa- 
rameters entering the model, nor to account for nonequi- 
librium exciton energy distribution and cooling due to 
phonon emissioniiii. Instead, we adopt a more phe- 
nomenological approach and assume that, as a result of 
cooling, the system can be described by an effective tem- 
perature, which leaves the densities n Ci h,x as the only im- 
portant hydrodynamical variables. The sources J c and Jh 
in Eq. describe the carrier photo-production, as well 
as the leakage current in the QW structure. 

In general, one can expect the electron-hole binding 
rate w and, to a lesser extent, the exciton recombination 
rate 7 to depend sensitively on the local exciton density 
n x . Of the several mechanisms that could lead to such a 
dependence at low temperatures close to exciton degen- 
eracy, perhaps the most important in the present context 
involve stimulated electron-hole binding processes medi- 
ated by phonons. These processes enhance the binding 
rate by a factor / = 1 + N^, where iV^ q denotes the oc- 
cupation of exciton states. In thermal equilibrium, and 
at low temperatures, one can ignore the reverse processes 
of exciton dissociation: The binding energy, carried away 
by phonons, is much larger than k^T. 

For a degenerate exciton gas with Ne=o > 1, the dom- 
inant process involves scattering into the ground state 
and the stimulated enhancement factor is expressed as 



no(T) 



n T = y ; 2 B . 2 



Here m x ~ 0.21 mo represents the exciton mass, and 
g denotes the degeneracy (in the indirect exciton sys- 
tem, the exchange interaction is extremely weak, and 
g = 4). Equivalently, when reparameterized through 
its dependence on temperature, u = Tq/T where Tq = 
(irh 2 /2<?to x /cb) n x is the degeneracy temperature. At 
T ~ Tq (equivalently n x ~ no), there is a crossover 
from classical to quantum Bose-Einstein statistics, and 
the stimulated enhancement factor / increases sharply. 

Qualitatively, the stimulated transition mechanism for 
hydrodynamic instability can be understood as follows: 
A local fluctuation in the exciton density leads to an in- 
crease in the stimulated electron-hole binding rate. The 
associated depletion of the local carrier concentration 
causes neighboring carriers to stream towards the point 
of fluctuation presenting a mechanism of positive feed- 
back. The wavelength, determined by the most unstable 
harmonic of the density, characterizes the length scale of 
spatial modulation in the nonuniform state. 



Before turning to the analysis of instability, it is use- 
ful to discuss intrinsic constraints on the dynamics i[TJl 
due to electric charge and particle number conservation. 
These are obtained by considering the linear combina- 
tions (c) — (h), (e) + (h) + 2(x) of the transport equa- 
tions In both cases, the nonlinear term drops out 
and one obtains linear equations 



with L, 



L c n c 



(h,x) 



+ L h n h 
dt-D, 



Je — Jh 



2jn x = J c + Jh 



(3) 
(4) 



Note that, since the ori- 



gin of the relations l^il 111 is routed in conservation laws, 
they are robust and insensitive to the exact form of the 
electron-hole binding term. 



B. The two-dimensional problem 

The simplest to analyze is instability of a steady state 
with spatially uniform carrier sources and density dis- 
tribution, which can be achieved in the presence of spa- 
tially extended photoexcitation, described by constant 
J c (r) = Jh(r) = J in Eqs.JTJ. In a uniform steady state, 
ignoring the dependence of the radiative recombination 
rate 7 on exciton density, we have 



= J/7, n e ,h = (J/w{n x )) 1/2 . 



(•5) 



The stability of the system can be assessed by linearizing 
Eqs. (JIJ about the uniform solution JSJ) with a harmonic 
modulation <5n e ,h,x oc e At e lk r . Using (l.'il 111 and writing 
L e 5n e = LhSrih — —L x Sn x — 7<5n x , one can express Sn^ x 
in terms of Sn Q and obtain 



L e (A,k)+ 7 ^ 1 



L c (X,k) \ = juLejX, k) 
L h (A,k)y L x (A,k)+ 7 



(0) 



where u = dlnw/dhin x is evaluated at the steady state 
©, and L c(h:X) (A,k) = A + /J e(hjX) k 2 . Solving Eq. ®, 
one obtains the growth rate dispersion A(k). The stabil- 
ity criterion Re A < is first violated at the parameter 
values such that the sign of A(k) reverses for some k. The 
condition A(k) = gives 



(k£ x ) 2 



(k4 



1 + (k4) 2 



(7) 



where £ x = \J D x /j denotes the exciton diffusion length 
and r = D^D^ 1 + DZ )n x /n e . Eq. J2J) has solutions if 



dhxw 
din n x 



> u r 



1 



,1/2 



(8) 



as illustrated in Fig. ^ At u — u c , we obtain the most 
unstable wavenumber fc* = r 1 / 4 jl x selected by competi- 
tion of the stimulated binding and diffusion processes. 

The binding rate w(n x ) builds up near degeneracy due 
to the growth of stimulated processes, leading to insta- 
bility at temperatures approaching Tbec- (For a weakly 
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interacting Bose gas, the transport coefficients, and thus 
the constant r, are practicaily insensitive to the degree 
of exciton degeneracy at T > Trec-) 
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FIG. 1: Graphical solution of Eq. that selects the most 
unstable wavelength. Inset: The 3-fold symmetric star of 
wavevectors describing the modulation near the instability 
threshold and the corresponding triangular pattern of exciton 
density variation. 

To what extent are these results insensitive to the ori- 
gin of the nonlinearity in the binding rate? If enhanced 
by intraband Auger processes, which transfer the binding 
energy released in exciton formation to other excitons, 
one expects the binding rate w to scale linearly with lo- 
cal exciton density, viz. w(n x ) = wo(l+n x /no), where no 
denotes some constant involving a ratio of the two-body 
and three-body cross-section of the electron and hole in 
the presence of excitons. Crucially, in this case, the left 
hand side of Eq. (JSJ is bounded by unity, while the right 
hand side is in excess. Therefore, at least over the param- 
eter range considered here, one can infer that a simple 
linear scaling of the binding rate with density does not 
lead to instability. Indeed, the instability may be used to 
discriminate against certain mechanisms in the kinetics 
of exciton formation. 

Turning to the discussion of the spatial pattern result- 
ing from the instability, we note that the wavevector se- 
lection determines its modulus, but not direction. At 
threshold u = u c all modes with |k| = /c* become unsta- 
ble simultaneously. The resulting 2d density distribution 
can be found by considering the effect of mixing differ- 
ent harmonics due to higher order terms in expanded 
in <5n e , x about the uniform state. Since these equations 
contain quadratic terms, the favored combination of har- 
monics is a 3-fold symmetric star 

k« =fc,(cos(f i + 0),sin(f j + 0)), (9) 

j = 1,2,3, with the phase parameter 9 describing the 
degeneracy with respect to 2d rotations. This leads to 
a density distribution Sn cx Y]j e ±kj ' r with maxima ar- 
ranged in a triangular lattice. 



On symmetry grounds, since the triangular lattice pat- 
tern is stabilized by quadratic terms, the mean field anal- 
ysis predicts that the transition to the modulated state in 
this case is abrupt, of a type I kind. Indeed, the triangu- 
lar lattice geometry arises in various 2d pattern selection 
problems, from Benard convection cellsi^ to the mixed 
state of type II superconductors 19 . 

C. Instability of a ID electron-hole interface 

The application of these ideas to the Id modulation 
seen in exciton rings^ 3 - was explored recently in Ref. 20 . 
The analysis starts with first determining the profile of 
the uniform distribution. The rings mark the inter- 
face between regions populated by electrons and holes 
at which they bind to form excitons. The steady state 
is maintained by a constant flux of carriers arriving at 
the interface. The parameter regime which is both rel- 
evant and simple to analyze is that of long exciton life- 
time 7 _1 where the diffusion length £ x exceeds the range 
of the electron and hole profile overlap. In this case, 
approximating the source of excitons by a straight line 
c5(x), where c is the total carrier flux and x is the coor- 
dinate normal to the interface, the exciton density profile 
is given by (c4/2L> x )e~l x l/^. Accordingly, one can seek 
the electron and hole profile treating w(n x ) as constant 
and restoring its dependence on n x later when turning 
to the instability. The profiles can be inferred from two 
coupled nonlinear diffusion equations 

<9 t n c (h) - D c{h) d 2 x n c(h) = -wn c n h , (10) 

with the boundary condition: D^d x n^\± a<1 = ±c8(±x). 
From Eq. @ one obtains D c n c — D^n^ = cx which allows 
the elimination of n^. Applying the rescaling n e (h) = 
c ^S'e(h)/-Ce(h)i where £ = (DeD-^/wc) 1 ^ 3 , one obtains 

dj g e = g c {g c - x), x = x/i. (11) 

From the rescaling one can infer that the electron and 
hole profiles overlap in a range of width I ~ c -1 / 3 while 

g c (\x\ > t) = x6(x) + O (Jir 1 /^- 2 ^ 372 / 3 ) . 

The solution of this problem, obtained numerically, is 
displayed in Fig. [3 

Although Eqs. IjlOfl are nonlinear, their diffusive char- 
acter does not straightforwardly admit a spatial instabil- 
ity: A fluctuation in the position of the interface initiates 
an increased electron-hole flux which, in time, restores 
the uniform distribution. However, if one restores the de- 
pendence of the binding rate on exciton density, the same 
mechanism of positive feedback which characterized the 
instability in the uniform system becomes active. 

To explore the instability, one may again expand lin- 
early in fluctuations around the spatially uniform solu- 
tion, g c (x, y) = g c (x) + Sg c (x) e tky (similarly g h and # x ), 
where y is the coordinate along the interface and g e (x) 
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FIG. 2: Numerical solution of the (normalized) bound state 
wavefunction "4>(x) representing fluctuation in the electron 
density, together with the effective potential U c g = g B + <?h — 
a o ffcgh of the Schrodinger-like equation H'lti . 

denotes the uniform profile obtained from Eq. . With 
£ x 3> £, the exciton density remains roughly uniform over 
the electron- hole interface. Denoting this value by n x (0), 
in the vicinity of the interface, one may again develop 
the linear expansion u>[n x ] ~ u>[n x (0)] (1 + uSn x /n x (0)) 
where, as before, u = dlnw/dlnn x . Noting that Eq. © 
enforces the steady state condition Sg c = Sg^, a lineariza- 
tion of Eqs.Q obtains the Schrodinger-like equation 

+ {£kf + g c + g h ] Sg c + ^-rg e 9h Sg x = 0, (12) 

where the particle density Fourier components 
<55x(o,h)(g) = / e~ iqx Sg x{c ^ h) (x)dx obey the relation 



where a\ = dx g e {x)g\ l {x)tp 2 {x) ~ 0.254^, and 
02(Q) = => / dxdx' g e (x)g h (x)i>(x)e-W x - x ^^{x'). 

* J -oo 

With £ x > £, it will follow that Q£ < £ x , and the latter 
takes the constant value a 2 — 0.461£ 2 independent of k. 
Finally, minimizing Su c with respect to k, one finds that 
the instability occurs with a wavevector 

fc c 4^ (a a 2 £ x /2£ 3 ) 1/3 (14) 

implying a shift of u c by Su c /u c ~ (^/^ x ) 4 ^ 3 - As a re- 
sult, one can infer that the spatial modulation wavelength 
A c ~ £ 1 / 3 ^ x ^ 3 is typically larger than the electron-hole 
overlap £, but smaller than £ x . Finally, an expansion 
of the nonlinear equations to higher order in fluctuations 
shows that, below the transition (i.e. for u > u c ), the am- 
plitude of the Fourier harmonic k c grows as (u — Uc) 1 ^ 2 . 

Once the instability is strongly developed (or when 
£ x < £) the above linear stability analysis fails to produce 
quantitatively accurate results. Instead, one can perform 
numerical search for a steady state of the full nonlinear 
problem. Such an analysis, reported in Ref^i, indeed 
demonstrates a modulational instability developing in a 
soft way, with a square root dependence of the mod- 
ulation amplitude on the deviation from critical point 
(see Fig. 3 in Ref."^). The characteristic instability wave- 
length obtained numerically is consistent with that pre- 
dicted by Eq. ltT4l . 



- -*9M (i - cyTt^) ■ (13) 

with Q 2 = (k£ x ) 2 + 1. Now, since the product g e gh is 
strongly peaked around the interface, the typical con- 
tribution from the last term in (|12|l arises from Fourier 
elements q£ ~ 1. Then, with £ x ^ £, the second contribu- 
tion to 5g x (q) can be treated as a small perturbation on 
the first and, to leading order, neglected, i.e. 5g x ~ — Sg c . 
In this approximation, the most unstable mode occurs at 
k = 0. Qualitatively, an increase in u will trigger an in- 
stability of the k = mode at a critical value u c when 
the linear equation first admits a non-zero solution for 
Sg c . At the critical point, the corresponding fluctuation 
in the electron density then acquires the profile of the 
(normalized) zero energy eigenstate ip(x). Numerically 
one finds that the critical point for the instability oc- 
curs when u c /g x {0) — (2£/£ x )u c = aa ~ 6.516, while the 
corresponding solution ip(x) is shown in Fig. [2 

While the approximation above identifies an instabil- 
ity, a perturbative analysis of the fc-dependent correc- 
tions implied by l|13(l reveals that the most unstable mode 
is spatially modulated. To the leading order of perturba- 
tion theory, an estimate of the shift of u c obtains 

u c a ai \ £$Q J 



III. LONG-RANGE ATTRACTION 
MECHANISM 

Here we discuss the hydrodynamical forces between ex- 
citons that arise due to polarization charge relaxation 
triggered by exciton recombination. For the indirect ex- 
citons in 2d coupled quantum wells we find a long-range 
attraction that falls off very slowly as a function of dis- 
tance. The strength of this interaction is controlled by 
the specifics of the model, such as carrier mobility in the 
wells and adjacent doped regions. Even if week, this in- 
teraction may overcome the repulsive electric dipole in- 
teraction at a certain length scale, thereby making the 
uniform state unstable with respect to density modula- 
tion formation. 

The idea of the 'plasmon wind' effect is as follows. 
Each exciton polarizes charges in the doped regions above 
and below the quantum wells. In the upper and lower 
doped region, the polarization charges are of opposite 
signs, forming a dipole that partially cancels the field of 
exciton dipole. While this partial screening effect is small 
due to a large distance to the doped region, and can be 
safely ignored in the static picture, it acquires another 
dimension and leads to a new effect if one takes into ac- 
count that excitons recombine at a constant rate. 

When an exciton recombines, its dipole disappears and 
the polarization charge, suddenly left to its own, starts 
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spreading in 2d. The spreading is controlled by charge 
continuity and electro/hydrodynamics of the doped re- 
gion. Thus an exciton cloud in which excitons disappear 
at a constant rate will represent a steady source of 2d 
plasmons spreading radially away from it. The polariza- 
tion charge in this situation is described by a 2d Poisson 
equation with a source proportional to exciton density 
and, therefore, the force on other excitons far away will 
be of a long-range 1/r form. The sign of this force is at- 
traction, since the polarization charge dipole is opposite 
to that of exciton dipole. 

The significance of this force is that, while being weak 
compared to the direct exciton repulsion, it falls off slowly 
and comes to dominate the physics at sufficiently large 
distances. Being attractive, it favors exciton clumping 
on large length scales. The situation can be understood 
by analogy with gravitation forces which are very weak, 
but long range, and thus are irrelevant on small length 
scales, while controlling the structures forming on the 
astronomical length scales. 

To estimate the effect of dynamical screening and plas- 
mon wind on the interaction between indirect excitons, 
we consider a simplified model. The excitons will be 
treated as electric dipoles positioned in the xy plane and 
oriented along the z axis. The doped regions will be de- 
scribed as 2d conducting planes with conductivities <Ji,2, 
at a distance ai,2 from the xy plane and parallel to it. 

The dynamical screening problem for indirect excitons, 
described by dipoles d with density n x (r, t), is formulated 
in terms of the charge and current densities in the doped 
regions: 

d t Pl,2 + Vjl >2 = Sl,2 , jl,2 = -0"1,2 V^l,2 (15) 

where s\ = — X(pi — P2)/2, S2 = — s\ describes leakage 
through the structure. The potentials v?i,2 are related to 
the charge density in the two planes as 

<fH = ^2u ij (k)p j + <p[ a) (16) 
i 

with the exciton dipole potential in the planes 

z = a\ 1 —a2- Using the Coulomb interaction 2d Fourier 
components 

t'-'/fts^"' (17 » 

we obtain 

U l=] (k) = 2ir/k , V m (k) = (2n/k)e- k ^ +a ^ (18) 

Similarly, <^ 0) = 2Trde~ kai , tp { 2 ° ] = -2ixde~ ka ^ is found 
for the potential in the planes z = Oi, —0,2 of a dipole d 
situated at the origin. 

To determine the force on a remote exciton, f = 
— V£/efij U e g — —dE z , we compute the electric field z 
component E z (k) — — 2ir(e~ kai p\ — e~ ka2 P2) using p\fi 



obtained from Eqs. i|15|) . l|16[l . l n the simplest symmetric 
situation, when o\^. — & and = we obtain 

U cff (k, w) = 2^ 2 ( - 2fc ^ 12 . - + k) n(k, w) 
\2k z Uii + (A — iuo)/<j J 

(19) 

where Un = U\\ — f/ 12 and n(k,uj) is exciton density. 
The second term corresponds to the direct exciton dipole 
interaction, while the first term describes the effect of 
polarization charge. 

In order to extract the static interaction, as well as 
the part that depends on the exciton recombination, we 
expand Ea. (|19|l to the first order in u> and obtain 

U cB (k,Lj) = (D^k) +iwD 2 (*))nx(*,w) (20) 

with 

™ = ^( wu*Z + V. + k ) (21) 

(22) 

At large distances, taking the small k limit of the expres- 
sions l(T5jl. i~ x -C k -C a" 1 , we obtain 

D 2 (k) = ^, D 1 (k)=(3 (23) 

where a = d 2 /(8aa 2 ), (3 — Trd 2 /a, and I = (Snaa/X) 1 ^ 2 is 
the characteristic length of plasmon decay due to vertical 
leakage across the quantum wells. The term D2, which 
has a long-range form — ln(|r — r'\/£) in position rep- 
resentation, describes the plasmon wind effect. The term 
D\ describes the static inter-exciton interaction screened 
by charges in the doped regions. Altogether, we obtain 
an expression 

UcS (r) = ^ x (r) - ^ / In (^^) M?')d 2 r' 

"" ^ ' ' (24) 

valid at the length scales from a to I. In passing from 
Ea. (|20|) to Ea. H24l) . we replaced n x by — 7n x , with 7 the 
recombination rate. (We note that a similar plasmon 
shakeup accompanies exciton formation by electron and 
hole binding, leading to qualitatively similar inter-exciton 
attraction.) 

To see how the long-range attractive part of U e g gen- 
erates instability, we consider exciton dynamics 

d t n x + Vj x = -jn x + J (25) 
j x = -D\7n x - n x n x \7U cS (26) 

where p x is exciton mobility. In a simple, although some- 
what artificial situation when the excitons are generated 
uniformly over the area, the homogeneous state is char- 
acterized by constant density n x = J/"/. Linearizing 
Ea. (|25|) about this state, with Sn x cx e"'e lkr , and set- 
ting lo < 0, we obtain the instability criterion: 

(D + p x n x f3)k + 7 < (27) 
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This condition is simplest to analyze in the absence of 
leakage, t = oo. Then for instability one must have 

a/i x n K > 1 (28) 

Replacing conductivity in a — d 2 /(8aa 2 ) by e 2 /^"-d, 
where \id is the mobility and rid is the carrier density 
in the doped region, we obtain the criterion 

d 2 fJ, x n x > 8e 2 a 2 (i d nd (29) 

In a coupled quantum well structure with well separation 
w the dipole moment is given by d 2 = e 2 w 2 /ecaAs where 
£GaAs is the dielectric constant. 

The criterion 1)29(1 shows that the instability can be 
reached by increasing exciton number n x over a thresh- 
old value n XiC = 8eGaAs(a/w) 2 (/W^ x )n d . In tne geome- 
try 1 ^, the ratio a/w ~ 10. With the factor 8ecaAs — 10 2 , 
even for the in-plane mobility /i x orders of magnitude 
higher than the mobility in the doped region, the critical 
density n KiC is quite high. 

While this difficulty might be overcome by employ- 
ing a stronger long-range force, e.g. due to phonons, a 
more serious issue with this mechanism arises from the 
consideration of temperature dependence. In the experi- 
ment!^, the instability occurs at temperature below few 
Kelvin. To explain this, the above model would have to 
account for strong temperature dependence of transport 
properties such as that taking place in the low conduc- 
tivity regime of electron hopping transport. In constrast, 
the kinetic effects discussed above lead in a natural way 



to a low temperature phase transition which occurs as 
a precursor of exciton degeneracy without any specific 
assumptions about the nature of transport. 



IV. CONCLUSION 

In summary, we considered a scenario when the real- 
ization of quantum degeneracy in a cold electron/hole- 
exciton system is signalled by the development of a spa- 
tial density modulation. Although our discussion is mo- 
tivated by the photoexcited CQW system in which elec- 
trons and holes are spatially separated, the mechanism is 
quite general applying also to geometries where the elec- 
tron and hole sources J e ,h have a spatially independent 
profile. The instability mechanism appears to depend 
sensitively on there being a strongly nonlinear depen- 
dence of the electron-hole binding rate on the exciton 
density pointing to the importance of stimulated scatter- 
ing. Although on a qualitative level this model is consis- 
tent with experimental observations of a low temperature 
transition, before considering it as the only explanation 
more experiments are needed in order to provide a quan- 
titative assessment of the instability mechanism. 
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